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Abstract 

We  cite  techniques  that  immediately  yield  theoretically  efficient  algorithms  for  computing 
solutions  to  eigenfunction  problems  for  self-adjoint  linear  differential  operators  of  any  finite  order 
acting  on  real- valued  functions  on  the  interval  [—1, 1];  these  algorithms  incur  computational  costs 
that  are  nearly  optimally  small,  to  within  factors  that  are  constant  multiples  of  small  powers 
of  the  logarithm  of  the  problem  size.  However,  the  factors  in  the  computational  costs  are  in 
general  too  large  for  practical  applications;  the  algorithms  would  appear  to  require  careful, 
detailed  optimizations  for  each  particular  application  in  order  to  be  useful  in  practice. 

For  any  positive  real  number  e,  positive  integer  m,  and  matrix  A,  we  say  that  A  is  semiseparable 
to  (absolute)  precision  s,  with  a  semiseparability  rank  of  at  most  m,  to  mean  that  any  contiguous 
rectangular  block  of  A  that  does  not  contain  any  entry  from  the  diagonal  of  A  has  at  most  m 
singular  values  greater  than  e. 

Observation  1  The  adaptive  versions  of  the  discretizations  described  in  [2],  [5],  and  [7]  recast 
eigenfunction  problems  for  self-adjoint  linear  differential  operators  of  any  finite  order  acting  on 
real- valued  functions  on  the  interval  [—1,1]  as  eigenvector  problems  for  self-adjoint  semiseparable 
matrices  (provided  that  all  boundary  conditions  are  linear  whenever  any  boundary  conditions  are 
prescribed) . 

Observation  2  For  any  positive  real  number  e,  positive  integers  m  and  N,  and  N  x  N  matrix 
A  that  is  self-adjoint,  has  N  distinct  eigenvalues,  and  is  semiseparable  to  precision  e,  with  a 
semiseparability  rank  of  at  most  m,  the  techniques  introduced  in  [4] ,  [3] ,  and  [1]  immediately  yield 
well-conditioned  algorithms  for  computing 

1.  the  N  eigenvalues  of  A,  at  a  cost  of  0{m  ■  N  -In  N)  operations  and  memory  storage  elements, 

2.  the  N  entries  of  any  of  the  eigenvectors  of  A,  at  a  cost  of  0{m  ■  N), 

3.  the  N  entries  of  any  linear  combination  of  the  normalized  eigenvectors  of  A  from  the  coeffi¬ 
cients  of  the  eigenvectors  in  the  linear  combination,  at  a  cost  of  0{m  ■  N  ■  InN),  and 

4.  the  N  coefficients  of  the  eigenvectors  in  any  linear  combination  of  the  normalized  eigenvectors 
of  A  from  the  entries  of  the  linear  combination,  at  a  cost  of  0{m  ■  N  ■  In  A); 
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these  costs  pertain  to  computations  performed  to  precision  e  and  with  A  itself  stored  efficiently,  at 
a  cost  of  0{m?  •  •  In  N). 


As  is  implicit  in  [1],  combining  Observations  1  and  2  immediately  yields  theoretically  efficient 
algorithms  for  computing  solutions  to  eigenfunction  problems  for  self-adjoint  linear  differential  op¬ 
erators  of  any  finite  order;  these  algorithms  incur  computational  costs  that  are  nearly  optimally 
small,  to  within  factors  that  are  constant  multiples  of  small  powers  of  the  logarithm  of  the  prob¬ 
lem  size.  However,  the  factors  in  the  computational  costs  are  in  general  too  large  for  practical 
applications.  Techniques  such  as  those  described  in  [6]  reduce  the  factors  for  specific  applications. 
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